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Developing Processing Techniques for Skylab Data 
Monthly Progress Report, March 1975 


The following report serves as the twenty-fifth monthly progress 
report for EREP Investigation 456 M which is entitled "Developing 
Processing Techniques for Skylab Data". The financial report for 
this contract (NAS9-13280) is being submitted under separate cover. 

The purpose of this investigation is to test information extraction 
techniques for SKYLAB S-192 data and compare with results obtained in 
applying these techniques to LANDSAT and aircraft scanner data. 

The SKYLAB S-192 data set being studied under this contract is the 
same data set being studied by a group from Michigan State University under 
the direction of Dr. Lester V. Manderscheid. Inasmuch as ERIM has sub- 
contracted with MSU to perform S-192 data analysis under contract NAS9-13332, 
a monthly report similar to this reporting the common data preparation is 
also being issued under that contract. 

During the reporting period we continued to emphasize the data pre- 
paration aspects of the task in an effort to finish entirely this part 
of the work. The decision to perform most of the processing on the non- 
scan line straightened, or conic, format data (which ia discussed below) 
caused this process to be drawn out further ..han had been anticipated. 

Jobs performed during the month included extending field location to conic 
scan line and point coordinates, the marking and digitization of field 
location points from the second U-2 acquired imagery, investigations into 
SDO to SDO misregistration in the conic and straightened data and measure- 
ment of pixel size for the S-192 data for the Michigan EREP test site data. 

To begin with, we were concerned over the effects of misregistration 
of SDOs on processing of S-192 data. A report [1] issued concerning S-192 
sensor evaluation called out 4 SDOs which were not perfectly registered. 
Examining the conic data, we found a sizeable water body where three scan 
lines made the transition from land to water at the 6atne point. Signals 
from these scan lines were averaged and the resulting data normalized to 
the boundary value and plotted. The results, shewn in Figure 1, showed the 
same misregistration problems as reported in the reference. We attempted 
a similar analysis on the straightened data, however we could not locate 
lake areas where the boundary occurred at the same point op several con- 
secutive scan lines. An analysis was made using just individual scan lines 
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FIGURE 1. MISREGISTRATION OF CONIC FORMAT S-I92 DATA 

(The symbols ) indicate the relative projection on the ground 
for all SDO's for one resolution element.) 
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but this failed because the noise in the data made the individual scan 
line traces too irregular to analyze. Ve could reach no conclusion 
regarding the misregistration of the straightened data. 

It was therefore decided to pursue a more analytical technique to 
measure the misregistration of the data. An algorithm was developed 
which calculates the cross correlation function for a pair of SDOs for 
many different simulated degrees of misregistration between the SDOs. 

The actual misregistration between the SDOs is that for which the cross 
correlation function is a maximum. The algorithm Is detailed in Appendix 1 
to this report. We plan during the coming month to implement the algorithm 
by programming it for our computer, and then using the program to measure 
the misregistration for both the conic and straightened data sets. 

However, it was felt that further processing should not be suspended 
pending the outcome of the above program. Since we wanted to continue 
with the processing effort and since we were convinced of the existence 
of misregistration, it was decided to continue the processing effort on 
the conic data. We feel that we could correct the conic data for mis- 
registration since the algorithm to correct for misregistration is simple 
for the conic data and we felt that we had a good estimate of the mis- 
registration of the conic data. We would further mention here that we 
have experienced very little problem in working with the conic data 
especially in regard to obtaining line and point numbers of particular 
pixels. Since the test area is located near the middle of the data swath, 
we have found that ordinary graymaps are only slightly distorted and art 
eminently useable for our needs. 

Having decided to process the conic data it was next necessary to 
convert the previously obtained straightened data line and point numbers 
for the fields in the ground information area to conic line and point 
numbers. This was done by using the inverse of the scan line straightening 
transformation equations as given in the EREP Users Handbook, coupled with 
regression techniques to accurately calculate the constants in the equations . 
The equations we used are: 
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CONIC POINT 
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P - [STRAIGHT POINT - 517.8-0.5] 

N ■ 1239 Points/Conic Scan Line 
0 - 116.25° Field of Scan 

A & B are constants estimated from regression techniques. 


Similarly, for scan lines: 


CONIC LINE = C + D • STRAIGHT LINE 


XI 


» E.R COS I (CONIC POINT * 2 - 2 ~ K)6 


with 


R * Radius of the scan circle projected on the Earth 
R % 608 pixels 


and C, D, and E are constants estimated from regression techniques. To 
perform the regression, 18 points were located on both conic and straightened 
graymaps. The regression fit was very good and further, all 5 coefficients 
seemed to be sensible, a reflection of the physical reality. 

With the field coordinates converted, the ground information was merged 
with the conic data. Graymaps of two conic data channels and the ground 
information channels were overlayed for comparison and the conversion was 
deemed very satisfactory. 

Finally, we ran the data through a deskewing program to reduce the 
misregistration in the data. With reference to Figure 1, we took the even 
numbered, high sample rate SDOs, alonf. with all the low sample rate SDOs. 

Thus there were three channels out of registration with the rest, two of 
these by a full pixel. The registration algorithm shifted these two channels 
over by one full pixel. The other SDO was registered by estimating the signal 
expected midway between the two samples. Initially the estimator used was a 
simple average of the two adjacent samples. 




WILLOW BUN LAtfOftATOftlC* TMt UNlVCfffttTV OF MICHIGAN 



101900-5 2-L 

Page 6 

During the reporting period we completed the marking of fields and 
other points of interest on the second U-2 acquired photograph, and 
digitised all these points as we had done for the first photograph. 

Again following the procedure established while working with the first 
photo, regression analysis was performed to calculate parameters to convert 
from photo coordinates to scan line straightened coordinates. Then the 
digitized points were converted. The next step is to convert, again using 
regression techniques to estimate parameters, to conic line and point 
numbers and to then merge the ground information with the data. We plan 
to do this during the next month. 

Other work during the period included the following three items: 

1) We received screening film of all SDO's (except for 15 and 16 
which are redundant) and examined it thoroughly. We expect this to be 
a useful tool throughout the analysis of the data. 

2) A brief analysis of actual pixel (not reselm) size was conducted. 
Fairs of pixels in lakes were located on straightened data graymaps which 
were either on the same scan line., several hundred pixels apart, or 
located at the same scan point number several hundred scan lines apart. 

Points corresponding to the pixels were also located on USGS maps of 
Southern Michigan. Distances were accurately measured on the USGS maps 
and on the graymaps with the result that the pixels were measured to be 

69 meters wide along scan direction and 72 meters in the along track 
direction. Calculations based c»n geometrical considerations using only 
the angle of the scan cone and the altitude at the time of data acquisition 
yielded measures of 70 x 70 meters. The EREP Users Handbook calls out the 
pixel size to be 72 x 72 meters. The differences are not felt to be serious. 

3) Crop acreage was totaled lor each of the 90 sections in the ground 
information area and grand totals were also calculated. This showed the 
predominant ground covers and percentage of total area as: 


Corn 

30% 

Trees 

17% 

Gras -ts 
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Stubbles 

9% 

Soil 

7% 
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During the next reporting period* in addition to those items already 
mentioned above, we intend to extract spectral signatures and perform 
analyses of the signatures, including calculation of optimum bands. With 
a set of signatures we will then he able to begin work on applying the 
mixtures classifier to SKYLAB data. We will also be able to evaluate 
signature extension techniques as applied to SKYLAB data. Further, we 
also intend to study the effects of misregistration on classification 
results for scan line straightened data where misregistration may not 
be correctable. 


Reference : 

[1] "ERIM Contributions to the S-192 Sensor Performance Evaluation", 
John G, Braithwaite and Peter F. Lambeck, ERIM 1Q2800-51-F, 
January 1975. 
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APPENDIX 1 


The following is a procedure for determining the amount of mis- 
registration between two correlated data channels. By reconstructing 
the continuous waveform over a lengthy interval in both channels, the 
cross-correlation function of the two waveforms can be determined. Let 
f(t) and g(t) denote the reconstructed waveforms in the two channels 
over the interval [A,C]. The cross-correlation function r(t ) is defined 
as ° 


C 

r(t Q ) - | f(t)g(t + t Q ) dt 

A 


The amount of misregistration between the two channels can be estimated 
as the value of the parameter t 0 which maximises the cross-correlation. 
The continuous waveforms can be reconstructed from the sample values by 
making assumptions which allow the use of Shannon's sampling theorem. 

The sampled data is converted into continuous form to allow the misregis- 
tration to be estimated to within a fraction of a pixel rather than in 
whole pixel increments. Hie length of the interval [A,C] must be long in 
comparison to the range of the parameter values t Q . This condition is 
required to minimize the effect of inaccuracies which will occur near the 
endpoints of the interval. 

Shannon's sampling theorem indicates that a continuous signal y(t), 
bandlimited to B (radians/sec) , can be exactly reconstructed from samples 
taken with a sampling interval t = -n/B. The sampling rate is equal to 
twice the highest frequency component contained in the signal. The 
original signal y(t) can be expressed in terms of the sample values y(mi) 
as 

y<t> - t B r(m> - T > 

m=-« 




Assume that the two continuous data channels f(t) and g(t) are band- 
limited to B and that the sampling interval t is equal to w/B. Let the 
sample values of these two waveforms over the interval [A,C] be denoted 
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«■ f(ki) and g(it) i, k * 1, ... , N. The croas-correlatlon r(t ) can 
be expressed in terns of the samples as ° 




N 

Y *<kT)g<iT) 
i,k-l 


l 


sin B(t - kt) sin B(t + t rt - It) 

B(t - ki) B(t + t - ii) 

o 


dt 


Using a variation of Parseval's Theorem, the integral can be evaluated 
by extending the limits of integration to positive and negative infinity, 
and r(t Q ) can be expressed as 


Jig sin B(kt - it + t ) 

r(t Q ) - Bit ^ f(kt)g(it) °“ 

i,k»l 


B(kt - it + t^) 


or, since Bt c it 


N 


— ^ sin ir(k - i + t /t) 

r(t 0 ) » B. £ t(kT)g(10 , (k - 1+ t 

i,k=l 0 


Thi? 8 relationship can be expressed in terms of a fraction of a sampling 
interval (or fraction of a pixel) by defining a variable A = t Q /t. Then 

r(4) - B. £ f(kt)g(l,) ~ * + *> ■ 

i,k=l 

Neglecting the constant factor Bti and expressing f(kT> and g(it) as f fc 
and g^, respectively, the function d(A) must be evaluated, where 


N 

iw- £ 

i,k-l 


f sin tr(k - 1 ■f A) 
k 8 i r(k - i + A) 


2p 
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which can be simplified as 

v l) r 

d(A) • I 

J— (N-l) L 1 

For large N, the variable j need not extend over the entire range because 
of the insignificant contribution of the high magnitude terms. To reduce 
the effects of noise, the function d(A) should be determined for several 
scan lines and averaged. 


-k-J 


f k *i 


n(A - j) 


